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ABSTRACT 

We perform Halo Occupation Distribution (HOD) modeling to interpret small-scale and 
intermediate-scale clustering of 35,000 luminous early- type galaxies and their cross-correlation with a 
reference imaging sample of normal galaxies in the Sloan Digital Sky Survey. The modeling results 
show that most of these luminous red galaxies (LRGs) are central galaxies residing in massive halos of 
typical mass M ~ a few times 10^^ to l{y^^h~^MQ^ while a few percent of them have to be satellites 
within halos in order to produce the strong auto-correlations exhibited on smaller scales. The mean 
luminosity Lc of central LRGs increases with the host halo mass, with a rough scaling relation of 
Lc oc M°'^. The halo mass required to host on average one satellite LRG above a luminosity thresh- 
old is found to be about 10 times higher than that required to host a central LRG above the same 
threshold. We find that in massive halos the distribution of L^, galaxies roughly follows that of the 
dark matter and their mean occupation number scales with halo mass as M^'^ . The HOD modeling 
results also allows for an intuitive understanding of the scale-dependent luminosity dependence of 
the cross-correlation between LRGs and galaxies. Constraints on the LRG HOD provide tests to 
models of formation and evolution of massive galaxies, and they are also useful for cosmological pa- 
rameter investigations. In one of the appendices, we provide LRG HOD parameters with dependence 
on cosmology inferred from modeling the two-point auto-correlation functions of LRGs. 
Subject headings: cosmology: observations — galaxies: halos — galaxies: statistics — galaxies: clus- 
ters: general — galaxies: elliptical and lenticular, cD — galaxies: evolution 



1. INTRODUCTION 

The clustering of galaxies depends 
erties, such as m o rphology (e.g. , 
Zwickv et all [1961 lOa-vis fc Gelled 
1980; P ostman k Gelleii 11984 1: iGuzzo et al 
Willmer. da Costa fc Pellegrinil [I998I : 
20021: iGoto et all I2003D . luminosity (e.g. 
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iMadgwicket al.l 120031 ). 




are typically the central regi ons of group s and clust ers. 
The LRG redshift sample Hisenstein efall 120011) o f 
the Sloan Digital Sky Survey (SDSS; ,York et all 12000) 
provides an enormous data set with which to measure 
the clustering of LRGs. In this paper, we perform 
theoretical modeling of auto-clustering of LRGs and 
cross-clustering of LRGs with other types of galaxies 
to understand the origin of their clustering properties, 
to learn how they are distributed among massive dark 
matter halos, and to aid the study of formation and 
evolution of massive galaxies. 

The theoretical understanding of galaxy cluster- 
ing has been greatly enhanced through the frame- 
work of the halo occupati on dis tribution ( HOD , 

see, e.g., Jing, Mo, & Borner 1998; Ma 



Fry ""200 



Galaxy clustering thus pro- 
vides important clues to the physics of galaxy formation. 
Often found to reside in galaxy groups and clusters, 
luminous red galaxies (LRGs) constitute the bright end 
of the galaxy luminosity function. Clustering of LRGs 
encodes information about their environments, which 
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Peacock & Smith 2000; Seljak 2000; Scoccimarro et al. 
20011: iBerhnd fc Weinberel l2002) and the closely related 



approach of the conditional lumin osity function (CLE, 
I Yang. Mo. fc van den BoschI l2003f l. The HOD formal- 
ism describes the bias relation between galaxies and 
matter at the level of individual virialized dark matter 
halos, whose distribution and properties can be readily 
predicted by numerical simulations or analytic models 
given a cosmological model. The key ingredients of 
this formalism are the probability distribution P{N\M) 
that a halo of mass M contains N galaxies of a given 
type and spatial and velocity distributions of galaxies 
within halos. In the CLE approach, the dependence 
of P{N\M) on galaxy luminosity is implicitly derived 
by inferring the conditional luminosity distribution 
of galaxies as a function of halo mass M . Given the 
HOD/CLE and the cosmology, any statistics of galaxy 
clustering on any scales can be predicted, therefore the 
HOD/CLE method provides a complete description of 
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the bias relation between galaxies and dark matter. 
HOD/CLF modeling has been applied to interpret 
galaxy clustering measurements in several g alaxy sur 
veys (e.g., Ijing & Borner 1998; Jing. Bor ner. fc Sut 
200l [Bo llock. Wechsler. fc Somervillel 1200 

Moustakas fc So mcrvillc 200^ : [ van den Bosch et al.l 



2001" 
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Yan. Ma dgwick. fc White..2003l: lZhens |l200l 'Yan g et al 
2005; Zeh avi et al.ll2005hHLee et aHl2 006: Haman a et al 
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Conrov et al. 
fc Zehavi 2007: 



200 &: I White et al 
Bla ke et all 120081 : 



Wake et al.ll2008( ). HOD/CLF analysis recasts galaxy 
clustering measurements into a form that is more 
physically informative and conduci ve for testing galaxy 
formation theories (see, e.g., iB erlind fc Weinberg! 
2002; Berlind et al. 2 0031: Ivan den B osch et al. 2003; 
Zheng et a l. 2005). ,Zehavi et al.l ( 2005bl) present cluster- 
ing measurements and HOD modeling for galaxies in the 
SDSS main galaxy spectroscopic sample ptrauss et al.l 
l2002f ). The HOD modeling results of the luminosity 
and color dependence of galaxy clustering are found to 
be in good qualitative agreement with predictions from 
galaxy formation models (Zheng et al. 2005). 

In this paper, we apply HOD modeling to inter- 
pret clustering in the SDSS LRG spectroscopic sample, 
which uses color and magnitude cuts to effectively se- 
lect LRGs out to redshift z ~ 0.45 (jEisenstein et al.l 
l2001h . The large volume probed by the LRG sample has 
led to the detection of the baryon acoustic peak in the 
two-point co r relation function (EisensteirLgt al. 20P5lj) . 
IZehavi et al.l (|2005aD report the measurements of two- 
point correlation functions of 35,000 LRGs on scales of 
0.3-40/i-iMpc. They find that LRGs are highly clus- 
tered (correlation length ~10ft-~^Mpc) and that more lu- 
minous LRGs are more clustered. Clear deviations from 
a power law are seen in the correlatio n functi ons, with 
a dip at ^2h~^y{^c. lEisenstein et al.l ()2005aD measure 
the two-point cross-correlation between 32,000 spectro- 
scopic LRGs and 16 million galaxies in the SDSS imaging 
sample. Since these reference galaxies have luminosities 
aroun d i», the characteristic luminosity of the Schechte"? 
(119761) luminosity func t ion, th ey are denoted as L, galax- 
iesl [Eisenstein et al.l ()2005aD find a strong luminosity 
dependence of the LRG-L* cross-clustering amplitude. 
The form of the luminosity dependence is itself depen- 
dent on scale, with more variation in the clustering am- 
plitude on small scales. Understanding all these auto- 
and cross-clustering features is one of the goals of our 
HOD modeling. Since LRGs trace massive halos, the 
cross-correlation between LRGs and I/* galaxies also al- 
lows us to study the HOD of galaxies in these massive 
halos. 

The structure of this paper is as follows. In § 2, we 
briefly describe the SDSS LRG and galaxy samples 
we use and our modeling method. In § 3, we describe our 
HOD parameterization for LRG samples and L, sam- 
ples. In § 4, we perform HOD modeling for LRG two- 
point auto-correlation functions and the two-point cross- 
correlation functions between LRGs and galaxies. We 
show how LRGs occupy dark matter halos and how 
galaxies occupy massive halos. Based on the modeling 
results, we interpret the luminosity dependence of the 
cross-clustering. Finally, we summarize and discuss our 
results in § 5. We also include three appendices in the 



paper. In Appendix [XI we investigate which HOD pa- 
rameter plays the major role in the departures of the 
galaxy two-point correlation function from a pure power 
law. In Appendix [Bl we present HOD modeling results 
for two luminosity threshold LRG samples for different 
cosmological models. In Appendix [C] we have a brief 
discussion on the mass function of the most massive ha- 
los and the fluctuation of the number of massive halos in 
the volume probed by the SDSS LRG samples. 

2. SAMPLES AND METHOD 

The LRG two-point auto-c orrelation functions i n the 
SDSS have been measured by IZehavi et al.l (|2005aD with 
spectroscopic samples. The luminosity cuts of the two 
LRG samples at z ~ 0.3 we model in this paper are 



-23.2 <Mg < 



-21.2 and 



-23.2 < Mg < 



-21.8, where 



Mg is the restframe g-band absolute magnitude at z = 
0.3 computed from the observed r-band magnitude with 
k and passive evolution corrections (jEisenstein et al.l 
[2005al) . Since galaxies with Mg < —23.2 are extremely 
rare, these can be regarded essentially as luminosity- 
threshold samples. For brevity, we call them the Mg < 
—21.2 sample and the Mg < —21.8 sample, or simply the 
faint sample and the bright sample. The comoving num- 
ber densities of the two samples are 9.73 x 10~^/i^Mpc~^ 
and 2.40 x 10~^h'^Mpc~^ , respectively. Because of the 
fiber collision effect, the smallest scale can be probed 
by the spectro s copic LRG samples is about 0.4/i~^Mpc. 
iMasiedi et al.l (|2006l ) extend the measurements of the 
two-point correlation function for the Mg < —21.2 sam- 
ple down to a scale of ~0.015/i~^Mpc by cross-correlating 
the spectroscopic sample with the imaging sample to 
avoid the fiber collision effect. They find that the real- 
space two-point correlation function of the Mg < —21.2 
sample roughly follow an profile from ~100/i~^Mpc 
down to ~0.01/i~^Mpc. In this paper, we focus on 
the measurements from the spectroscopic samples in 
IZehavi et al.l (|2005all and limit ourselves to modeling the 
clustering on scales above 0.3ft.~^Mpc. 

The two-point cross-correlation functions between 
spectroscopic LRG samples and the imaging galaxy 
sam ple in the SDSS at z ^ 0.3 are measured by 
lEisenstein et al.l ([2005a). The galaxy sample we 
model is the one defined by the luminosity range M* — 0.5 
to M* + 1.0, w h ich is called the M* + 1.0 sample in 
lEisenstein et al.l ()2005af ). For the LRG sample in the 
cross-clustering, we adopt the luminosity-bin sample de- 
fined by -21.7 < Mg < -21.2. 
We essent i ally f ollow the method presented in 
(l2005bl) and adopt the improvements of 
(l2005f ) for theoretical modeling of the two- 



Zehavi et al 
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point auto- and cross-correlation functions in the HOD 
framework. In this paper, more general HOD parame- 
terizations are used, as described in the following sec- 
ti ons. While the two -point auto-correlation functions 
in IZehavi et al.l (|2005af l are p rojected along the red shift 
direction, similar to those in IZehavi et al.l (j2005b( ). the 
LRG-L* two-point cross- c orrelat ion functions A mea- 
sured in 'Eisenstei n et al.l (|2005a[ ) are volume-averaged 
real-space cross-correlation functions ■ The average 
is weighted by a spherical window function. 
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where the window function is of the form 



^2 / ,.2 

W^(r;a) = -exp U — 



(2) 



The effective volume V for the window function is 
3(27ra^)'^/^. In our modeling, we perform the same vol- 
ume average to compute the predicted cross-correlation 
function. 

When performing fits to the two-point auto- and 
cross-correlation functions, we calculate values of us- 
ing the full error covariance matrice s, inferred thro ugh 
the iackknife meth od (see Eisenstein "eran[2005al and 
IZehavi et aLll2005bl for details). 

Throughout the paper, we adopt the spatially flat 
"concordance" cosmological model with the matter den- 
sity parameter flm = 0.3 and baryon density parameter 
rib = 0.047. We assume adiabatic Gaussian primordial 
density fluctuations with a power-law index of the spec- 
trum Us = 1. The r.m.s. matter density fluctuation 
in spheres of radius 8ft.~^Mpc linearly extrapolated to 
z = is assumed to be as = 0.8. The Hubble constant 
we use is h — 0.7 in unit of 100 km Mpc~^. Ap- 
pendix B presents the dependence of the derived HOD 
parameters on erg and Qm for a five-parameter HOD 
model. The mean redshift of the above LRG and L» 
samples is around 0.3 and our model calculations take 
this into account, i.e., all halo properties are calculated 
at z = 0.3 for the adopted cosmology. In general, the 
comoving u nit is adopted fo r distan ce, but to be com- 
patible with ' Eisenstein et al.l (|2005aD . we also use proper 
units with explicit mention when discussing the cross- 
clustering. Finally, we assume h — 1 when quoting mag- 
nitudes throughout the paper. 

3. HOD PARAMETERIZATION 

3.1. The Mean Occupation Function of LRGs 

The two LRG samples for which the two-point auto- 
correlation functions are modeled are nearly luminosity- 
threshold samples. The LRG sam ples defined in 
the cro ss-correlation measurements in lEisenstein et al.l 
(j2005ar ) are luminosity-bin samples. We concentrate on 
the cross-correlation between L^, galaxies and —21.7 < 
Mg < —21.2 LRGs, together with the two-point auto- 
correlation functions of the two luminosity-threshold 
samples with Mg < -21.2 and Mg < -21.8. The 
three LRG samples with different luminosity cuts pro- 
vide leverage to constrain the relation between galaxy 
luminosity and halo mass. We therefore first parameter- 
ize the LRG HOD in a way that includes the luminosity 
information. 

For HOD parameterization, it has been found to be 
useful to separate contribution s from central and s atellite 
galaxies (jKravtsov et all 12004 : IZheng et al.ll2005[ ). Cen- 
tral galaxy luminosity is correlated with the host halo 
mass. Bas ed on predictions of galaxy formation models 
() Zheng et a l. 2005), we assume that at a fixed halo mass 
the central galaxy luminosity follows a log-normal distri- 
bution and the mean luminosity of this distribution has 
a power-law form Lc = Ls{M /MsY , with p independent 
of halo mass in the range that LRGs probe. In terms of 
absolute magnitudes. 



where Mgs is the mean luminosity in halos with a pivot 
mass Ms and is simply set to be —19.8. We note that in 
this paper M with a subscript "(;" stands for the g-band 
absolute magnitude, while that without the subscript "g" 
is used for halo mass. With a constant standard devia- 
tion CTMg in magnitude, the luminosity distribution of 
central galaxies in halos of mass M is then 



rf(iVeen(M)) 
dMa 



1 



27rf7 M, 



■ exp 



[Mg - Mg^f 



(4) 



LRG samples are not defined purely through luminosity, 
since there is a color cut in the LRG selection. There- 
fore, for an LRG sample with Mg < Mg^thros, the mean 
occupation function of central galaxies is calculated as 



(7Vec„(M)) 



M„ 



dM, 



d{N,UM)) 

dMn 



fiMg), (5) 



Mg, = Mgs - 2.5plog{M/Ms), 



(3) 



where f{Mg) is the LRG selection function at a given 
luminosity. Based on the LRG luminosity function and 
the total galaxy luminosity function, we approximate the 
selection function as /(Mg) = -2(Mg + 20.45)/3 in the 
range -21.95 < Mg < -20.45, for Mg > -20.45, and 1 
for Mg < —21.95. The HOD of a luminosity-bin galaxy 
sample is just the difference of those of two luminosity- 
threshold samples. 

The theoretically predicted mean occupation function 
of satellite galaxies fo r a luminosity-t hreshold sa mple is 
close to a power law (jKravtsov et al.l , 2004: Zhen g et al.l 
I2005D . Here, we adopt a more flexible parameteriza- 
tion, which can give us a better idea on the constrain- 
ing power of the two-point clustering on the HOD. For 
each luminosity-threshold LRG sample, we parameterize 
the HOD of satellite LRGs through the mean occupation 
numbers at five mass scales, and the the mean occupa- 
tion function is assumed to be a cubic spline curve pass - 
ing through the five points (see lZheng fc Weinbergl2007f l. 
Linear extrapolations are used outside of the ranges of 
the five points. In addition, we impose the central galaxy 
cutoff profile (iVccn(M)) to (Nsi^tiM)). The five mass 
scales are chosen to be log{M/h-^MQ) = 14.10 -I- 0.35i 
(i — 0, 1, 2, 3, 4). The mean occupation function of satel- 
lites in a luminosity-bin sample is obtained by interpo- 
lating the values of \og{Nsat{M)) for the two luminosity- 
threshold samples (linearly in log L) . The satellite occu- 
pation number aroun d the mean is assumed to follow the 
Poiss on distribution (jKravtsov et al.l 120041 : iZheng et al.l 
l2005l) . 

The parameterization of the halo occupation of LRGs 
we present here is a hybrid of the CLF framework 
(I Yang. Mo. fc van de n Bosch 200.^ and the HOD frame- 
work (Ber lind fc Weinberg! I2002f l — the distribution of 
central LRGs is put in a CLF form, while that of the 
satellite LRGs is in the usual HOD form. This combina- 
tion allows us to infer information on how luminosity of 
central galaxies changes with halo mass and at the same 
time to keep a reasonable number of free parameters. 

With the above parameterized HOD, we can eas- 
ily form the mean occupation functions for the two 
luminosity-threshold LRG samples and the luminosity- 
bin sample given their luminosity cuts. The thirteen free 
parameters are the pivot mass scale Ms and the slope p 
in the luminosity-mass relation, the scatter ctm of the 
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luminosity distribution of central LRGs, and the total of 
ten spline points of {Nsat(M)) for the two luminosity- 
threshold samples. In Appendix [Bl we provide HOD 
modeling results for the two luminosity threshold sam- 
ples and their cosmology dependence based on a simpler 
five-parameter HOD model. 

3.2. The Mean Occupation Function of Galaxies 

The sample in lEisenstein et al.l ()2005af ) is at a 
mean redshift z ~ 0.3. To parameterize the HOD of 
this Lj, sample, we make use of the modeling result in 
IZehavi et al.l (|2005bD at z ~ 0. The L* galaxy sample is a 
luminosity-bin sample M*-0.5 < Mr < M* + 1.0, which 
can be regarded as the difference between two luminosity- 
threshold sa mples, Mr < M * + l.O and Mr < M* - 0.5. 
As shown in iKravtsov et al.| (|2004f) . at low redshifts, for 
threshold samples of a fixed (comoving) galaxy num- 
ber density, the shape of the mean occupation func- 
tion as a function of redshift approximately remains un- 
changed. We take the advantage of this presumed prop- 
erty and construct the mean occupation function for 
galaxies at z ~ 0.3 based on the HODs of z ^ 
galaxies. We first find out the number densities of the 
two luminosity-threshold samples with Mr < M* -f 1.0 
and Mr < M* - 0.5 at z ~ 0.3. To do this, we 
ad opt the evolut i on m odel of the luminosity function 
in iBlanton et al.l (|2003f) set to z = 0.3 and obtain the 
two number densities as ni = 1.1 x 10~^/i^Mpc~'^ and 
712 = 1.3 X 10~^/i^Mpc~^, respectively. We find the lu- 
minosity thresholds Li and L2 for z '--^ galaxies that 
match the two number densities, 

n(L > Li;z - 0) = ni (6) 

and 

n{L > L2; z 0) = n2. (7) 

We then use the results presented in'Zehavi et al.' (2005b') 
to infer the HODs at z ~ for luminosity-threshold sam- 
ples with Li and L2 as the thresholds. Finally, we take 
the difference of these two HODs to infer the mean oc- 
cupation function for galaxies in the Li < L < L2 lumi- 
nosity bin, 

{N{Li <L< L2]M)) = {N{L > Li;M))-{N{L > L2]M)). 

We take the shape of this mean occupation function as 
that for the galaxies at z ~ 0.3. The form of this 
mean occupation function is a sum of a square window 
for central galaxies and a power law with an index 1.10 
for satellites. With the shape fixed, we allow the mean 
occupation function to shift in mass scales. Furthermore, 
we add freedom to the high mass slope, as described be- 
low, because we have better data to constrain it. 

As shown later, the small-scale two-point cross- 
correlation between galaxies and LRGs comes from 
satellite L^, galaxies paired with LRGs. As our goal 
here is to investigate the occupation in high mass 
halos, it is necessary to introduce additional degrees of 
freedom in the mean occupation function of the satel- 
lite galaxies in massive halos. We only adopt the 
z ~ 0.3 mean occupation function for satellites, as 
constructed above, up to lO^^/i~^M0, and introduce an- 
other power law to represent the mean occupation func- 
tion of (satellite) L* galaxies in halos more massive than 



IQ^^h^^MQ. We assume this occupation distribution fol- 
lows the Poisson distribution with the mean in the form 
of Nt{M /Mt)°'*^ . Here is a pivot mass scale fixed to 
be 2.5 X \Q^'^h~^MQ^ Nt is the mean occupation num- 
ber of -L* galaxies in halos of this pivot mass, and 
denotes the slope at the high mass end. The scale of the 
pivot mass is chosen to minimize the correlation between 
Nt and oih- In order to have more flexibility, we do not 
impose continuity at lQ^^h~^MQ. 

At first glance, our construction of the i» HOD seems 
complex, involving theoretical priors and introducing an 
artificial break in the satellite occupation function. A 
more general HOD model for galaxies could be ob- 
tained by parameterizing the HODs of two luminos- 
ity threshold samples and taking their difference (e.g., 
[Tinker et al. 2007). However, we only have a small num- 
ber of data points from the LRG-L, cross-correlation 
functions to constrain the HOD. The constraints on 
the low mass part of the i» HOD (especially the mean 
occupation function for central galaxies) come from the 
the number density of the sample and the large scale 
bias of the cross-correlation. The shape of the low mass 
part of a general HOD can not be well constrained by 
either the number density or the large scale bias, given 
that the halo bias factor in this mass range only increases 
slowly. The constraints to the high mass part of the L* 
HOD come from the small-scale cross-correlation with 
the LRGs. Our procedure of the i» HOD parameteri- 
zation reduces the degrees of freedom in the low mass 
part of the HOD and allows some flexibility in the high 
mass part of the HOD which is our focus. Although 
our parameterization is a restricted one, it is suitable 
for our purpose of understanding clustering properties in 
L,-LRG cross-correlation. To summarize, we allow the 
low mass end of the mean occupation function (a square 
window plus the low-mass power law) to have an overall 
horizontal shift with the shape fixed. The main role of 
shifting the low mass part of the L* HOD is to match the 
number density of the L, sample and the large scale bias 
of the cross-correlation, while the high mass part of the 
occupation function of i» galaxies that we are interested 
in determines the small-scale cross-correlation with the 
LRGs. 

The distribution profile of satellite galaxies in- 
side halos is assu med to follow the Navarro - Frenk - White 
(NFW) profile (|Navarro. Frenk. fc Whit j [l995L II996L 
|1997|) with a concentration parameter c = co{M/M^\Y , 
where M^\ ^ 1.5 x IO^^/i^^Mq is the nonlinear ma ss at 
z ~ 0.3 and (3 is fixed to -0.13 (i Bullock et al]|2001h . As 
a whole, there are four free parameters in the L* HOD, 
the overall horizontal shift of the low mass part, the high 
mass amplitude Nt and slope ah of the mean occupa- 
tion function, and the halo concentration parameter cq 
extrapolated to the nonlinear mass scale. To compute 
the one-halo term of LRG-L* cross-correlation, we make 
the simplifying assumption that the occupation number 
of LRGs and that of L» galaxies inside the same halo 
are uncorrelated. That is, (iV(LRG; M)A^(L,; M)) = 
(7V(LRG;A/))(iV(i,; Af)). We discuss the effect of this 
assumption on our modeling results in § 4. 

4. MODELING RESULTS 

We model simultaneously the two-point auto- 
correlation functions of the two luminosity-threshold 
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Fig. 1. — Projected two-point auto-correlation functions and best-fit HODs for the two luminosity-threshold LRG samples. Panel (a): 
The measured two-point correlation functions (data points and error bars) and the HOD model fits (solid curves). The two dashed curves 
for each sample show the envelope of predictions from models with Ax^ < 4. The predicted one-halo and two-halo terms (dot-dashed 
curves) are also shown for the sample with the lower luminosity threshold. Panel (6): The mean occupation functions (solid curves) of 
LRGs from the best fits, with contributions from central (dotted) and satellite (dashed) LRGs. For each sample, the two sets of curves 
are the envelope from models with Ax'^ < 4 and those for all galaxies are shaded. Panel (c): The marginalized distribution of the satellite 
fraction in each sample with the central 68% distribution marked by the two dashed lines. Panel (d): The probability distribution of halo 
masses for the LRGs in each sample (solid lines), obtained from the occupation function shown in panel (6) weighted by the differential 
halo mass function. The dotted lines show the halo mass probability distribution from just the central galaxies. 



LRG samples {Mg < 



-21.2 and Mg < 



-21.8) and the 



two-point cross-correlation function between the and 
the luminosity-bin (-21.7 < Mg < -21.2) LRG sample. 
Altogether, there are seventeen free parameters, four for 
the HOD and thirteen for the LRG HOD. 

In addition to the sum of from the two auto- 
correlation functions and the cross-correlation function, 
we also add the number densities of the two threshold 
LRG samples and the galaxies into the overall x^- 
That is, 



= (wi - Wi)'^Ci ^(wi - wl) 

+ (W2 - W;)'^C2 ^(W2 - W2) 



(A 



E 



{rii 



(9) 



where wi, W2, and A are the vectors of auto-correlation 
functions of the two LRG samples and the LRG-L* cross- 
correlation function, and [i = 1,2,3) are the three 
number densities. The observed values are denoted with 
a superscript *. The full covariance matrix is used for 
each correlation function and 10% fractional errors are 
assumed for each of the galaxy number densities. 

The three covariance matrices, Ci, C2, and Cx 
are estimated using jackknife resampling. For each of 
the t wo LRG samples, 104 jackknife subsamples are 
used (jZehavi et al.. .2005a) for estimating the covari- 
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ance among 10 data points. For the LRG-L* cross- 
correlation, the covariance among 6 d ata points is es- 
timate d with 50 jack k nife su bsamples (jEisenstein et alj 
l2005al ). iZehavi et alJ ()2005bn performed extensive tests 
with mock catalogs to access the reliability of the jack- 
knife error estimates in projected correlation functions 
over a similar range of separations. They used 100 
mock catalogs with the same geometry and angular 
completeness as the SDSS sample and similar cluster- 
ing properties, created using the PTHalos method of 
[Scoccimarro & Shcth (2002). Their tests showed that 
the jackknife method is a robust way to estimate the er- 
ror covariance matrix, especially for the relatively large 
volumes probed here. This holds as long as the number of 
jackknife realizations, n, is signific antly larger than the 
dimension of the data vector, p. IHartlap et al.l (j2007T ) 
discusses related issues, pointing out a potential bias in 
general model fitting which depends on the ratio of p 
to n. As a crude test of our error uncertainties, we in- 
corporate their proposed method to remove this bias by 
multiplying the inverse of the covariance matrix by a fac- 
tor of {n — p— 2)/{n— 1). In our case, this factor is about 
0.9. We perform this for the Mg < 21.2 LRG sample, re- 
sulting in a ~ 10% increase in the uncertainties of the 
HOD parameters, but no noticeable change in the best- 
fit values. We note that the results presented later in 
this paper do not include such a correction. 

We assume the different clustering measurements to 
be independent, ignoring possible correlations between 
statistical errors across correlation functions and num- 
ber densities of different samples. Strictly speaking, 
such statistical correlations are not zero. For exam- 
ple, the Mg < —21.2 LRG sample includes LRGs with 
Mg < —21.8, therefore the auto-correlation functions of 
the Mg < —21.2 and Mg < —21.8 samples are partially 
correlated. These correlations could be estimated with 
the jackknife technique, but they would be noisy. Ne- 
glecting such correlations might make the constraints on 
the HOD parameters somewhat tighter than they should 
be, which is a caveat for interpreting our results. 

We adopt a Markov Chain Monte Carlo technique 
(MCMC; e.g., IGilks. Richardson, fc Spiegelhalteri [199^ 
to explore the HOD parameter space. At each point of 
the chain, we take a random walk in the parameter space 
to generate a new set of HOD parameters. The step-size 
of the random walk for each parameter is drawn from 
a Gaussian distribution. The probability to accept the 
new set of HOD parameters depends on the difference 
between ^^'^ Xncw i^^^ values of for the new and 
old models): 1 for xLw < xIm and exp[-(xLw - Xou)] 
for Xncw > Xo\d- F^^^ priors in logarithmic space are 
adopted for the LRG satellite occupation numbers at the 
ten spline points and the pivot mass scale Ms and for Nt 
and mass scale shift of galaxies. Flat priors in linear 
space are used for other HOD parameters. 

The total number of data points to model is 29 (10+10 
auto-correlation measurements for the two LRG samples, 
6 LRG-L* cross-correlation measurements, and 3 num- 
ber densities) . The number densities are well reproduced 
by our model, with the median values 1.1-a, 0.4-cr, and 
0.1-cr away from the observed ones for the two LRG and 
the L* samples. Our model has 17 free parameters and 
therefore the number of degree of freedom is 12. We find 



that the best-fit model has =25. The probability for 
a x^ v alue higher than our best-fit x^ is 1.5%. (Includ- 
ing the'Hartla p et all ()2007f ) correction gives x^ ~ 22 and 
the probability increases to ~4%). The relatively-large 
X^ values might be partly caused by our neglecting the 
correlation between statistical errors across correlation 
functions and number densities of different samples. It 
also indicates that the accuracy in our analytical model 
of the two-point galaxy correlation function needs to be 
improved and that our HOD parameterization is not per- 
fect. With the above caveats in mind, we present our 
modeling results below. 

4.1. Constraints on the HOD of LRGs 
4.1.1. HOD for the Luminosity-threshold LRG Samples 

Figure [1] shows the fitting results for the two 
luminosity-threshold LRG samples. Figure (TJa shows 
the best-fit projected two-point correlation functions to- 
gether with the measurements. For the Mg < —21.2 sam- 
ple, the one-halo and two-halo components (dot-dashed 
curves) of the fit are also shown. In Figure [TJ), the mean 
occupation function for each sample is plotted, separat- 
ing into contributions from central (dotted) and satellite 
(dashed) galaxies. The range of mean occupation func- 
tions with Ax^ < 4 is denoted by the shaded region 
(Ax^ < 4 is chosen so that the envelopes are sampled 
by a large number of MCMC points and the range delin- 
eated by the envelopes can be clearly seen in the plot). 
As in IZehavi et al.l (|2005bf ). the amplitude of the high 
halo mass end of the mean occupation function is poorly 
constrained by two-point correlation functions. Since the 
abundance of halos drops exponentially at the high mass 
end, the two-halo pairs arc mostly contributed by lower 
mass halos. Although the number of one-halo pairs per 
halo rises roughly as {N{M))'^, the exponential drop of 
the number of high mass halos also makes most one-halo 
pairs come from lower mass halos. Overall, the two-point 
correlation function is dominated by signals from halos 
of lower mass where {N{M))^ a few, leading to the poor 
constraints on the HOD at the very high mass end. The 
constraint on the high mass end slope of {Nsa,t{M)) for 
the faint LRG sample appears to be relatively strong, 
and we discuss its possible implications in § [SJ 

The results show that the host halos of LRGs are mas- 
sive, above 10^'^ft.~^Mo, corresponding to large galaxy 
groups and clusters, which is consistent wit h stu dies of 
the environment of luminous galaxies (e.g., iLohl 120031 . 
We find that as the threshold luminosity increases by a 
factor of ~1.7, the mass scale of host halos shifts by a 
factor of ~ 2.3, which implies a correlation between the 
luminosity Lc of the central LRG and the mass M of the 
host halo in the form of Lc oc M°-^^. 

LRGs are often thought to be the central elliptical 
galaxies of galaxy groups of clusters. Our results clearly 
show that a fraction of the LRGs have to be satel- 
lite galaxies. The main constraint for this comes from 
the high amplitudes of the correlation function at small 
scales. Without satellite LRGs, the projected correla- 
tion function would flatten out toward small scales, sim- 
ilar in shape to the two-halo term shown in Figure [TJa. 
To explain the rising small-scale clustering, the one-halo 
term has to be introduced, hence we require the exis- 
tence of satellite LRGs. However, the fits imply that 
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Fig. 2. — Constraints on HOD parameters from the LRG auto-correlation functions and the LRG-L* cross-correlation function (see the 
text for details). Panel (a): The marginalized distribution of the parameter crji/^ , which is the width (in magnitude) of central galaxy 
luminosity distribution at a fixed halo mass. Panel (6): The marginalized distribution of the parameter p, which characterizes the relation 
between mean central galaxy luminosity and halo mass through Lc oc M^. In panels (a) and (6), the two dashed vertical lines indicate the 
central 68% of the distribution. Panel (c): Illustration of the effect of crMg and p on the mean occupation functions of LRGs. The thick 
dotted and dashed curves are for the Mg < —21.2 LRG sample, with CMg varied from 0.37 to 0.43 (i.e., a ztlcr change). The thin dashed 
and dotted curves are for the Alg < —21.8 LRG sample, with p varied from 0.46 to 0.51 (a itlcr change in p). Panel (d): The marginalized 
joint distribution of the concentration parameter cq (normalized to that at the z ~ 0.3 nonlinear mass scale) and the high mass slope 
of the mean occupation function for the L* galaxies. The contours show the 68% and 95% confidence levels for two parameters. 



only a small fraction of LRGs are satellites: 5.2-6.2% 



and 2.3-3.2% in the Mg < 



-21.2 and Mg < -21.8 sam- 



ples, respectively, as shown in Figure[T]c. The steeper rise 
of the small scale clustering makes the two-point correla- 
tion function of the brighter LRG sample deviate from a 
power law more prominently. In Appendix |Al we investi- 
gate the key ingredients that lead to the departure from 
a power law and provide insight as to why the departure 
becomes more clear for samples that are more luminous 
(as shown here) a nd for samples at higher redshifts (e.g., 
lOuchi et al.l[2005h . 

The mean occupation function tells us the mean num- 
ber of LRGs as a function of the halo mass. The result 
can be viewed differently by asking what the probability 



distribution of masses of halos hosting such LRGs would 
be. The probability is simply the product of the mean oc- 
cupation function and the differential halo mass function. 
We show such probability distributions from the MCMC 
run in Figure [1^. The distribution of Af^ < -21.2 LRGs 
peaks at ~ 4.5 x IO^^H-^Mq, and for Mg < -21.8 LRGs 
the peak shifts to higher mass, ~ lQ^'^h~^MQ. Both 
distributions span a large range in halo mass. The full 
width at half maximum (FWHM) for either distribu- 
tion is about AlogM = 0.8. Since only a small frac- 
tion of LRGs are satellites, the probability distribution 
is almost determined by the central LRGs. Only in ha- 
los more massive than a few times lO^'^'h'^Mp, does the 
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number of satellites in a halo become significant (i.e., 
greater than one on average), and around this mass scale 
a low amplitude shoulder in the distribution emerges. In 
FWHM sense, we find that Mg < -21.2 and Mg < -21.8 
LRGs reside in halos of mass 2-13x 10^^/i~^Mq and 4- 
25xlO^^/i"^Mo, respectively. 

From the mean occupation functions in Figure [TJ*, one 
also notices th at, compared with t he samples of lower 
luminosities in iZehavi et al.l (|2005bf ) , the low mass cut- 
off profiles for the LRG samples are better constrained. 
The relatively tight constraints come from the steepen- 
ing both in the halo mass function and in the mass de- 
pendence of the halo bias factor toward high halo mass. 
The former steepening makes the galaxy number den- 
sity more sensitive to the cutoff profile, while the lat- 
ter one increases the sensitivity of the large-scale galaxy 
bias factor to the cutoff profile. With our HOD param- 
eterization, the cutoff profile of the LRG mean occupa- 
tion function encodes information on the distribution of 
central galaxy luminosity in halos of fixed mass and on 
how the mean central galaxy luminosity scales with halo 
mass. Figure [2]a shows the marginalized distribution of 
the width (in magnitude) of the (log-normal) central 
galaxy luminosity distribution at a fixed halo mass. The 
clustering data require a scatter of 0.16 dex in the cen- 
tral galaxy luminosity in halos of a given mass. The thick 
dotted and dashed curves in Figure [2]c shows how aMg 
affects the cutoff profile in (iVcon(M)) of the Mg < -21.2 
sample by varying ctm by ztltr. We note that the con- 
straints on the cutoff profile (and therefore ) depend 
on cosmological parameters, especially erg, in the sense 
of a larger for a larger dg (see Appendix |B] for more 
details). 

The mean central galaxy luminosity scales with halo 
mass, which is characterized by the parameter p, Lc oc 
MP. LRG samples with different luminosity cuts allow 
us to constrain the parameter p. The marginalized dis- 
tribution of p is shown in Figure ^j. The thin dotted 
and dashed curves in Figure [2]c show the effect p on the 
mass scale shift of the Mg < —21.8 sample relative to 
the Mg < —21.2 sample. The value of p is around 0.48, 
which seems to be inconsistent with the value 0.66 es- 
timated from comparing the luminosity and halo mass 
scale of the two threshold LRG samples. The reason is 
simple — in our parameterization, the parameters p and 
o'Mg constrained here correspond to all central galaxies, 
not only central LRG galaxies that have a color selection 
criterion imposed. 

Our HOD parameterization here is rather flexible in 
the satellite HOD. Again, in Appendix |B1 we present 
modeling results with a five-parameter HOD model and 
show their dependence on cosmology. In § [5l we compare 
our HOD modeling results with those from other works 
and discuss a few issues related to the modeling. 

We note that, on large scales, the bestfit Wp curve in 
Figure [T]a appears to be lower than that from the five- 
parameter model (Fig. [8h) for the Alg < —21.2 sam- 
ple. The low mass cutoff profiles of {Nccn{M)) for the 
bright and the faint LRG samples are correlated in the 
parameterization adopted here, which means that the 
HOD for central galaxies is more restrictive than that in 
the five-parameter model. The LRG-L.^ cross-correlation 
also limits the range of LRG HODs. The bestfit solution 



is a compromise in matching the large-scale amplitudes of 
both the LRG auto-correlation functions and the LRG- 
cross-correlation function. The result implies that 
our HOD parameterization is not perfect and that there 
is room to improve it (e.g., by allowing the scatter in the 
central galaxy luminosity to vary with halo mass). 

4.1.2. Mass Scales of Host Halos of Central and Satellite 

LRGs 

Applying HOD modeling t o galaxy samples wit h dif- 
ferent luminosity thresholds, IZehavi et al.l ()2005bl ) find 
that there is a remarkable scaling relation between the 
characteristic minimum mass Afmin of the host halos 
and the mass scale Mi of halos that on average host 
one satellite galaxy (above the luminosity threshold) in 
addition to the central galaxy. Mi ~ 23Afniin. With 
a HOD pa rameterization close to what is used here, 
iZheng. Coi l. & Zehavi (2007) found the relation is more 
like Ml 18Afniin. Given the parameterization de- 
scribed in ^'S.l[ we define Mmin to be the halo mass at 
which the expected number of central galaxies above the 
luminosity threshold is 0.5, whether or not the galaxy 
satisfies the LRG color criteria. Theoretical studies of 
HODs of subhalos in dissipationless dark matter simula- 
tions (Krav tsov et al. 2004) and those predi cted by SPH 
and s emi-analytic galaxy formation models (|Zheng et al.l 
l2005f ) reveal a similar r elation with the sca ling factor 
around 20. As shown bv lBerhnd et al.l ()2003[ ). the large 
gap between Mi and Mj^in arises because in the low oc- 
cupation regime, a more massive halo tends to host a 
more massive central galaxy, rather than multiple smaller 
galaxies. Does the Mi-M^nin scaling relation extend to 
massive halos hosting LRGs? 

In Figure [3l we plot Afn im and Mi as a function of th e 
threshold luminosity from lZheng. Coil, fc ZehavH (|2007[ ). 
corrected to be consistent with erg = 0.8 adopted in this 
paper. The luminositv in, Zheng. G oil. & Zehavi (2 00?! ) is 
in z = 0.1 r-band (see IZehavi et a l. 2005b). For compar- 
ison, we convert the i^-corrected and passively evolved 
z = 0.3 g-band threshold luminosities of the two LRG 
samples to the z = 0.1 r-band ones by adopting an ap- 
parent color g — r = 0.4, and we obtain Mr < —21.6 and 
Mr < —22.2, respectively. We obtain the distribution of 
Minin and Mi for the two luminosity-threshold LRG sam- 
ples modeled in this paper by solving {Ncen{M)) — 0.5 
and {Nsat{M)) = 1 for each set of HOD parameters in 
the MCMC chain. The results are shown as the last 
two pairs of points. Roughly speaking, they seem to 
follow the previous trend. The mild discontinuity proba- 
bly reflects an imperfect magnitude conversion, which is 
not surprising as we are trying to account for filter dif- 
ference, JiT-correction, and stellar population evolution 
from z 0.3 to z ~ 0.15. A more interesting difference 
is that the scaling factors between Mi and Mj^in become 
10.5l^;5 (Mg < -21.2) and 7.2l°;^ {Mg < -21.8), re- 
spectively. For comparison, the mean Mi-to-Mmin ratio 
for the MAIN galaxy samples shown in Figure [3] is 18.3 
with a mean la uncertainty of ~3.6. This indicates that 
the scaling relation may break down in very massive ha- 
los. In fact, the scaling relation from the MAIN sample 
with the highest luminosity threshold {Mr < —21.5) in 
Figure [3] has already shown a trend of such a decrease, 
with an Mi-to-Mmin ratio of 13.2^2 ^. As a whole, these 
results show that the A/i-Mmin scaling factor decreases 
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Fig. 3. — Mass scales of the LRG HODs as a function of thresh- 
old luminosity. Shown are the relation between the character- 
istic minimum mass Mmin at which 50% of halos host central 
galaxies above the luminosity threshold and the mass Mi of ha- 
los that on average host one satellite galaxy, as a function of the 
threshold absolute ma g nitud e. The square points are taken from 
IZheng. Coil, fc Zehavil II2007I ) for the SDSS main galaxy sample 
(corrected to erg = 0.8). Open and filled circles are the results 
for the two LRG samples (note that the z = 0.3 g-band luminos- 
ity is converted to 2 = 0.1 r-band by adopting an apparent color 
of g — r = 0.4). Dotted lines indicate different relations between 
luminosity of central galaxy and mass of host halo. 

for luminous galaxies that reside in massive halos. We 
caution, however, that the definition of the Mi/M^in 
ratio becomes more parameterization dependent at high 
galaxy luminosities because the scatter between luminos- 
ity and halo mass is larger (i.e., the low mass cutoff of 
the HOD is less sharp). 

The change of the Mi-Af,„i„ scaling factor in massive 
halos can be understood from the competition between 
accretion and destruction processes. In general, massive 
halos accrete their satellites more recently than less mas- 
sive halos. While the rate of satellite accretion for low 
mass halos peaks at ~ 10 Gyr ago, cluster-sized halos 
const antly accrete satellites until recently ( Zentncr 'et al.l 
120051 ). As a consequence, there is less time for the orbit 
of satellites in a massive halo to decay through dynamical 
friction and for them to merge with the central galaxy 
to form a larger (brighter) central galaxy. In addition, 
the LRG samples we study are at redshift ^ 0.3, which 
makes the accretion even more dominant. Therefore, the 
decrease in the ratio of Mi and M^i^ could be a man- 
ifestation of the favor of accretion over destruction in 
massive halos. 

4.1.3. HOD for the Luminosity-Bin LRG Sample 

The cross-correlation between luminosity-bin LRG 
sample (—21.7 < Mg < —21.2) and the galaxies also 
leads to constraints on the HOD of the luminosity-bin 
LRG sample. The lower (solid) curves in the right panel 
of Figure |4] are the Ax^ < 4 envelope of the mean occu- 
pation function of the -21.7 < Mg < -21.2 LRGs. The 
bump (dotted) around ^ lO^^/i^^M© is the contribution 
from central galaxies and the dashed curves represent 
contributions from satellites. The result shows that most 



of the LRGs in this luminosity bin are central galaxies in 
halos of mass ^2-20x10^^ Mq, and a small fraction 
(~7%) of them are satellites in more massive halos. 

4.2. Constraints on the HOD of Galaxies 

iZehavi et all ()2005bD perform HOD modeling of the 
two-point auto-correlation functions of the MAIN galaxy 
sample and show that the two-point function can im- 
pose important constraints on the HOD of a sample of 
galaxies. In general, the mean occupation function is 
tightly constrained around {N{M)) ~ a few. It becomes 
loosely constrained towards higher halo masses because 
the two-point correlation function is less sensitive to the 
occupation distribution in these halos as a result of the 
steep drop of the halo mass function. Therefore, the 
analyses of the MAIN galaxy s ample with low luminos- 
ity thresholds in IZehavi et al.l ()2005bf ) cannot quite re- 
veal how these low luminosity galaxies reside in massive 
halos. Since LRGs automatically pick out the massive 
halos, the cross-correlation between MAIN sample galax- 
ies and LRGs provides us a nice way to study the halo 
occupation of MAIN sample galaxies in massive halos, 
enabling us to better constrain the HOD of galaxies 
in massive halos. 

Figure Wfi shows constraints on the concentration pa- 
rameter Co and the high mass slope ah of the mean 
occupation function for galaxies, marginalized over 
the other parameters. These two parameters are cor- 
related in a sense that a higher cq corresponds to a 
higher ah- Higher ah means that more galaxies 
reside in higher mass halos with lower concentrations 
and larger virial radius, so to maintain the small-scale 
(cross-)clustering the distribution of galaxies need to be 
more concentrated, i.e., a higher cq. The high mass 
slope of the mean occupation function of galaxies is 
1.49 ± 0.09 (l-cr range), steeper than the value 1.10 at 
lower mass. The concentration parameter cq represents 
the value extrapolated to M^i ~ 1.5 x lO^^/i~^M0 ac- 
cording to c = co(M/M„i)-° i3. Its 1 -a range is found to 
be 9.4 ± 1.7, which translates to 5.4 ±1.0 and 4.0 ±0.7 in 
halos of lO'^'^h-'^MQ and IO^^H-'^Mq, respectively, im- 
plying that the distribution of L* galaxies more or less 
follows that of the dark matter, at radii the data can 
probe (> 0.2/i~^Mpc comoving). The mean occupation 
function of galaxies in massive halos is shown as up- 
per (solid) curves in the right panel of Figure 21 These 
two curves are the envelope determined by Ax^ < 4. On 
average, about ten galaxies are expected to reside in 
a halo of mass ^ lO^'^h'^ Mq. 

Our fits show that the mean occupation function of 
galaxies at high halo masses becomes steeper than 
the slope 1.10 at intermediate masses (and lower red- 
shift). We have also run a model with fixed slope 1.10 
for the satellite mean occupation function in the whole 
mass range. The fits to the LRG auto-correlation func- 
tion and the LRG-L^^ cross-correlation function become 
much worse, with the overall increasing by ~50. The 
resulting best fit to the cross-correlation function is plot- 
ted in the left panels of Figured) The lower panel shows 
the fractional difference between the fits (thick curve for 
the flexible model and thin curve for the one with fixed 
slope) and the data. The fit from the model with fixed 
slope is almost on top of the one from the more flexible 
model, but the error bars in the measurement are small 
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Fig. 4. — Mean occupation functions and cross-correlation functions of L* galaxies and LRGs from the HOD modeling. Left panels: the 
predicted two-point cross-correlation function between LRGs and L* galaxies separated into contributions from central and satellite LRGs 
paired with L, galaxies within same halos and LRGs paired with L« galaxies from different halos. The scale a is in units of proper h~^Mpc 
to be consistent with that adopted in iEisenstein et al.l Il2005al ). The quantity <j>o = 2.267 X 10~^/i^Mpc~'' is the proper number density of 
L, galaxies at the mean redshift z ~ 0.3. The lower panel shows the fractional difference between model fits (thick for the flexible model 
and thin for the slope-fixed model) and data (see the text). Right panel: the mean occupation function for the L» galaxies in massive 
halos (top solid curves) and that for the luminosity- bin LRGs. (bottom curves). The square window at the low mass end represents the 
mean occupation function of central L« galaxies. The envelopes of the mean occupation functions are derived from models with Ax^ < 4. 
Dotted and dashed curves in the LRG mean occupation function are contributions from central and satellite galaxies, respectively. The 
dot-dashed curve shows the mean occupation function for the L* galaxies from a model with the slope fixed to be 1.10 in the whole mass 
range. 

and the change in is substantial (an increase of ~ 10 
in the contributed by the cross-correlation function). 
The corresponding high-mass end of the mean occupa- 
tion function of galaxies from this more restricted 
model is shown as the dot-dashed line in the right panel 
of FigurelH The single-slope model has more satellites 
in halos less massive than IO^'^H'^Mq, and the best-fit 
concentration parameter for galaxies is about a factor 
of two smaller; these two effects compensate each other to 
approximately maintain the amplitude of the small-scale 
cross-correlation. Unless the error bars in the measure- 
ments were underestimated (e.g., by a factor of two), the 
restricted model is highly disfavored by its much worse 
fits to the data (Ax^ 50). 

The cause and implication of the steep inferred high 
mass slope of the galaxy occupation function is not 
clear. It may be related to the selection of galax- 
ies — the sample used here is composed of galaxies 
in a bin of 1.5 magnitude around a redshift-dependent 
reference magnitude that i s supp osed to match L» at 
z ~ 0.3 ([Eisenstein et al.l I2005af) . Since the redshift 
range (0.16 < z < 0.44) is not narrow, the sample 
should be regarded as an effective sample, rather than 
a uniform sample. It may also be caused by the im- 
perfection in the analytic model of the two-point cross- 
correlation function. The small measurement errors in 
the two-point cross-correlation function may require a 
more accurate model than the one used in this paper, 
and allowing for uncertainty in the model accuracy would 
increase the allowed range of the high mass slope. Fi- 
nally, it may also be related to the assumption that 
the occupation numbers of LRGs and galaxies in the 
same halo are uncorrelated (§ 3.2). Since the effect of 
any correlation between the two occupation numbers be- 



comes smaller as halo mass increases (jSimon "erall[2009l) . 
an anti-correlation between the LRG and occupation 
numbers would lead to a decrease in the high mass slope 
of the galaxy to conserve the number of LRG-L* 
pairs. As a whole, we therefore have higher confidence 
in the value of (N) - 10 at M ~ lO"/i"^M0, where the 
fits from the two models cross, than we have in the slope 
ah itself. 

The modeling result also shows (see the left panel 
of Figure H]) that the LRG-L* cross-correlation func- 
tion is dominated by central LRGs paired with satel- 
lite L■^, galaxies on scales less than 0.5/i~^Mpc (co- 
moving), while above this scale the signal from satellite 
LRGs paired with satellite galaxies takes over un- 
til ^ 1.5ft.^^Mpc (comoving), where the two-halo pairs 
start to dominate. We show below that the variation 
with scales in the contributing components is the key to 
understanding the scale dependence of the luminosity- 
dependent cross-clustering. In the one-halo term, the 
signal from central LRGs mainly comes from halos of 
mass a few times lO^'^ to 10^'^h~^MQ, while the signal 
from satelhte LRGs is from halos of - 10^'^-^ h~'^MQ, 
where the LRG occupation number is a few. Although 
the occupation numbers of LRGs and galaxies keep 
rising toward higher halo mass, more massive halos are 
too rare to make a significant contribution to the cross- 
correlation signal. 

4.3. On the Scale Dependence of the 
Luminosity-Dependent Cross- Clustering 

For the LRG-L* cross-correlation, lEisenstein et al.l 
(|2005al) find a strong dependence on LRG luminosity. 
The clustering amplitude becomes higher for more lumi- 
nous LRGs and varies by a factor of up to four over a fac- 



HOD Modeling of LRG Clustering 



11 



4 

3 

< 

>o 2 
-e- 

1 



15 



< 10 
> 



"T 1 — I 1 1 — I — I — I — I 1 — I 1" 

a = 0.125h-iMpc J 
- • Data |J' 

Illlllllllll Model j/ ^ i 



(a) 

_i I I I L 



2 4 6 8 
L/U 



"T — I — I — I — I — I — r~| — I — I — I — I — I — r 

" a = 0.5h-'Mpc 



t t 



-L- -3, 



illlllf ' ' 



■(., 

_l I I I '-j-' -I i" " "I" T "i" i ' T I I L. 




10 
8 

< 6 
> 

o 

■G- 4 
2 


20 
15 

> 
^ 10 

5 



80 
60 

< 

>o 40 



"T — I — I — I — I — I — I — I — I — I — I — I — r 

a = 0.25h-'Mpc 

■■■■ l-halo cen-sat A 

_ ■■■■ l-halo sal-sat if 

2-halo 



-y 



"(b) 

J I I I I 'f -I- -I ■ -I- - t - -I- - I - -I- - 1 • • 



2 4 6 8 

L/L, 



~i — I — I — I — I — I — I — I — I — I — I — I — I — I — r 

" a=lh"'Mpc T 



i 



(d) , 

J I I I ■l--l-"i'T I I I I I I L. 



20 




2 4 6 
L/L. 



~i — I — I — I — I — I — — I — I — I — I — I — I — r 

a=4h"'Mpc 



I T-i*'*" 



00 



II i 



(f) 



I . I ■ .1. ■ L ■ J- -I I I— J— -J— -J— .-I— 

2 4 6 8 

L/U 



Fig. 5. — The luminosity dependence of the LRG— L* cross-correlation functions at different scales. The six panels correspond to (proper) 
scales a= 0.125 to 4h~^Mpc, respectively, as labeled in each panel. The quantity <j>o = 2.267 X lO^^/i^Mpc""^ is the proper number density 
of L* galaxies at the mean redshift z ^ 0.3, and V = 3(27ra^)'^''^ is the effective volume for the window function [see eq. J2J]. In each 
panel, the shaded regions are predictions from the modeling results. Since the model predictions are based on clustering information up to 
only ~ 4L*, the shaded regions below and above 4L* are model interpolations and extrapolations, respectively. The data points with error 
bars are the measurements in [Eisenstcin ot al. (2005a). Thick (thin) dotted curves are contributions from central (satellite) LRGs paired 
with L* galaxies within common halos (calculated from the best-fit HOD model), and dashed curves, which can be clearly seen in panels 
(d)—{f), represent the two-halo pair contribution. 



tor of four in LRG luminosity (see their Fig. 2). Further- 
more, the clustering amplitude increases more strongly 
with luminosity at smaller scales. We now show that 
these complex trends can be largely explained by the 
HOD results described above. We note that the HOD 
results are based on LRGs with luminos ity L < AL^,, and 
while t he LRG-L* cross-correlation in lEisenstein et al.l 
(|2005a[ ) is measured up to an LRG luminosity of 8L*, we 
focus on the L < 4L* results. We emphasize that we do 
not intend to explain the data points for L > AL^, for 
which we extrapolate our results. Our purpose here is to 
give a qualitative interpretation of the scale-dependent 
luminosity dependence of the cross-correlation between 
LRGs and galaxies. 
From our modeling results, the HOD for LRGs in nar- 



row luminosity bins can be readily constructed, similar to 
what we do for the —21.7 < Mg < —21.2 sample. We do 
not refit auto-correlation functions of different LRG sam- 
ples; rather, we apply the value of (Tj\/^ and the scaling 
Lc oc MP derived from the faint and bright luminosity- 
threshold samples and interpolate (or extrapolate) the 
satellite occupation function from these samples. Fig- 
ure El compares the predicted dependence of the cross- 
correlation on LRG luminosity to the observed one. The 
plotted quantity (j)oVA is the average excess number of 
galaxies around an LRG in an effective spherical vol- 
ume of V (see eq. [l]). We have already shown (Fig. [4]) 
that our HOD model reproduces the overall scale depen- 
dence of the cross-correlation accurately. Figure [5] shows 
that the model also captures the trend of stronger lu- 
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minosity dependence at smaller scales. Since the predic- 
tion is mostly based on modeling LRGs in the luminosity 
range of 2-4 , it becomes less accurate at higher himi- 
nosities, where it is essentiahy an extrapolation. 

Close inspection of the observed points shows that, at 
the smallest scales, the clustering amplitude rises steeply 
and steadily from ~ 2L* to ^ 8L*; while on larger scales, 
the luminosity dependence in the 2L^ - AL^ range is rel- 
atively flat, before steepening at higher luminosities. In 
addition, the overall luminosity dependence is weaker at 
larger scales. An intuitive understanding of these fea- 
tures can be built from our HOD modeling results. 

At very small scales, over all the LRG luminosity range, 
the cross-correlation signal is dominated by central LRG 
paired with satellite L^, galaxies. The thick dotted curves 
in Figure [5] show this component from the model. The 
slope of this component is estimated as follows. Approxi- 
mating the mean occupation function of central LRGs in 
a narrow luminosity bin as a Dirac-(5 function, the two- 
point cross-correlation function at a separation r between 
galaxies and LRGs is simply proportional to the pair 
number {N^{M)) f{r; M), where {N^.{M)) is the mean 
occupation function of L* galaxies in halos of mass M 
and /(r; M) is the fraction of galaxies located at a 
radius r from the central LRG in halos of mass M. The 
function / is just the distribution proflle of galaxies, 
/(r; M) = p{r; M)r'^ / J p{r; M)r'^dr. Using a power law 
to approximate the inner profile, p(r; M) (x [r/R^-^)^, we 

have /(r;M) oc R~t^^'^^ oc Af-(i+T/3), where i?vir is the 
virial radius of the halo. Since {N^{M)) oc M"'', we see 
that the cross-correlation amplitude A oc M"''~^^+'''/3^. 
Noting that the central luminosity L oc M^, the cross- 
correlation amplitude has the dependence on luminosity 
as A oc Inserting typical values of the 

model results, an = 1.49, 7 = — 1 (inner NFW profile), 
andp = 0.66 to the expression, we obtain A oc i^'^, close 
to the observational result, which is roughly A oc L^'^. 

On larger scales (but still within the regime of one-halo 
pair domination), the contribution from satellite LRGs 
paired with satellite galaxies starts to show up at the 
low luminosity end, as can be seen at scales 0.25 and 
0.5 /i~-'^Mpc in Figure [51 At a fixed scale, most satellite- 
satellite pairs come from halos of a narrow mass range. 
As the luminosity in consideration increases, the occupa- 
tion number of satellite LRGs at this halo mass decreases, 
leading to a decreasing contribution to the clustering am- 
plitude from satellite-satellite pairs. The opposing de- 
pendences of clustering contributions from central LRGs 
and satellite LRGs on luminosity flattens the overall lu- 
minosity dependence at the low luminosity end, a feature 
seen in the observed clustering. 

On much larger scales (e.g., 4ft,~^Mpc), the two-halo 
pairs dominate the cross-correlation between LRGs and 
L, galaxies, and the signal is proportional to the large- 
scale bias factor of LRGs. If the LRG bias factor is 
approximated by the halo bias factor, the luminosity- 
dependent clustering simply reflects the dependence of 
the halo bias factor on the halo mass. Around M ^ 
10^"* /i~^M0, the halo bias factor can be approximated 
by a power law with index ah ~ 0.3-0.4. This, together 
with the relation between central luminosity and halo 
mass L oc {p ~ 0.66), gives a luminosity dependence 
of the cross-clustering following roughly as L°-^, which 



agrees well with the observed trend. 

In modeling the LRG-L, cross-correlation, we make 
the assumption that there is no correlation between 
the occupation numbers of LRGs and galaxies in- 
side the same halo. Loosing this assumption would lead 
to changes in the contributions from one-halo central- 
satellite and satellite-satellite pairs. These changes 
would be at the level of fine details, and the above pic- 
ture of the interplay among the three components for 
interpreting the scale-dependent luminosity dependence 
of the cross-correlation would remain valid. 

5. SUMMARY AND DISCUSSION 

We have modeled the two-point auto-correlation func- 
tions of LRGs and the two-point cross-correlation func- 
tions between LRGs and galaxies in the SDSS, within 
the HOD framework, obtaining results on the mean rela- 
tion between central LRG luminosity and halo mass, the 
dispersion about this relation, the slope and amplitude 
of the satellite occupation function, and the abundance 
of L, galaxies in massive halos. 

The continuous rise toward small scales of two-point 
auto-correlation functions of LRGs implies that not all 
LRGs are the bright central galaxies in galaxy groups or 
clusters; a fraction of LRGs must be satellites to pro- 
duce small scale, one-halo pairs. However, the satellite 
fraction is small and decreases with the LRG luminosity, 
e.g., ~5-6% for Mg < -21.2 and -2-3% for Mg < -21.8 
based on the HOD modeling. The characteristic mini- 
mum host halo mass of central LRGs (Mmin, at which 
50% of halos host a galaxy above the luminosity thresh- 
old) is a few times 10^^h~^MQ and increases with LRG 
l uminosity. 

IZehavi et all (|2005bl : see also IZheng. Cod, fc Zehavil 
found a ratio Mi/M^i^ ~ 20 between the halo 
mass required to host a satellite above a luminosity 
threshold and the mass required to host a central galaxy 
above the same threshold. For these LRG samples, which 
populate higher mass halos and have a median redshift 
z — 0.3, we find a smaller ratio, Mi/Afniin ^ 10. A 
similar drop is seen for th e brightest sample {M^ < 
—22) in IZehavi et al.l (j2005bD . which has a mean redshift 
~ 0.16. The decrease of the scaling factor reflects the 
balance between acc retion and destruction of satellites 
(|Zent ner "erall[2005h — massive halos assemble more re- 
cently and their satellites have less time to merge with 
the central galaxy. The relatively higher redshift of the 
samples further strengthens this effect. 

The HOD of LRGs has been inferred using different 
methods in several recent investigations. To compare our 
results with others, one needs to pay attention to the 
differences in the sample definition, the underlying as- 
sumptions, and the modeling details. Mandel baum et al.l 
(2006) present a mass determination of host halos for 
two SDSS LRG samples based on galaxy-galaxy lensing 
measurements. The construction of their LRG samples 
is not identical to ours, but the halo masses determined 
from their two samples appear to be consistent with those 
of the two luminosity-threshold LRG samples we model. 
The galaxy lensing directly measures the halo masses, 
while our results come from the mass distribution of ha- 
los and the galaxy assignment required to reproduce the 
observed clustering. The agreement between the two re- 
sults is therefore encouraging. .Wake et al.. (,2008il model 
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the projected two-point correlation functions for LRGs 
in the 2dF-SDSS LRG and QSO survey with a three- 
parameter model. For the SDSS z = 0.21 LRG sam- 
ple, their inferred mea n occupation functi on is in general 
agreements with ours. iBlake et all jlOOSl) perform HOD 
modeling of the two-point angular correlation function 
of 0.4 < z < 0.7 SDSS LRGs with photometric redshifts. 
Their parameterization is a slight variation of our five- 
parameter model presented in Appendix |B] For sam- 
ples with similar number densities, their inferred halo 
mass scales, cutoff widths of central galaxy occupation 
function, and high mass slopes of satel lite occupation 
functi on are close to what we obtain. iKulkarni et al.l 
((200l constrain the HOD of SDSS LRGs with redshift- 
space two-point and three-point correlation functions, by 
comparing the measurements to those from mock cat- 
alogs generated through populating halos identified in 
iV-body simulations. They use a three-parameter HOD 
description and assume no velocity bias. They find that 
redshift-space three-point correlation functions favor a 
lower high mass slope (~ 1.4) for the satellite occupa- 
tion function. Since the z = outputs of erg — 0.9 
simulations are used in their modeling while the median 
redshift of LRGs is about 0.3, the effective erg in their 
modeling is about 1.05. From Appendix [Bl (eq. [B3]), 
we see that the high mass slope from our modeling is 
1.5 for s uch a high q-g , which is close to the value fa- 
vored by IKulkarni et aP (|2007f l. However, we note that 
they use quite a different halo definition than ours (with 
a much lower overdensity threshold), which complicates 
the comparison. A more detailed comparison of the LRG 
modeling re sults with different methods and sam ples can 
be found in lBrown. Zheng. White, et al. I (|2008D . 

The inferred high mass slopes of the LRG occupation 
functions tend to be substantially lar ger than unity, ei- 
ther from our results or others fe.g., IBlake et al.l 120081 : 
IKulkarni et"ani2007l : IWake et al.ll2008[ ). This appears to 
differ from observational inferences and th eoretical pre- 
dictions for low luminosity samples (e.g., iZe havi et al.| 
2005b; Zheng, Coil, & Zchavi 2007; K ravtsov e t al. 2004; 
Zheng et al.ll2005HConrov et al.l i2006). Is the steep slope 
of the inferred LRG occupation function a true feature or 
merely a result of modeling imperfections? First of all, 
the steep slope should not be a result of any restriction 
in our HOD paramcterizations. Our five-parameter HOD 
model (Appendix |B| introduces a cutoff in {Nsnt{M)) at 
the low mass end so that the connection between the 
high mass end slope and that at (iVsat(Af)) ~ a few is 
broken. For the fainter LRG sample, the constraint on 
the slope up to M ^ 2 x IO^^H'^Mq remains tight even 
with our more flexible parameterization [see Fig.[T]J6)]. In 
the LRG survey volume, there are not many halos that 
are more massive than 10^^h~^MQ. The true halo mass 
function in this volume can therefore deviate from the 
theoretical one used in the modehng. In Appendix [Cl we 
show that this fluctuation in the halo mass function does 
not seem to introduce any systematic biases in model fit- 
ting as it is already reflected in the covariance matrix of 
the data. To determine which features of the data drive 
the steep slope, we fit the faint LRG sample with the 
five-parameter model after excluding some data points 
and find that the two data points at ~ 1.7/i~^Mpc 
and 2.7/i~^Mpc play a large role — if they are excluded. 



the slope drops from 1.8 to 1.6 (for erg = 0.8). These 
scales are in the one-halo to two-halo transition region 
where the model is sensitive to the treatments of halo ex- 
clusion and scale-dependent halo bias. A more accurate 
scale-dependent halo bias with halos defined by spherical 
over-density (SO) could lead to a somewhat lower value 
of the high r nass end slope ( J . L. T inker, private com- 
munication). iReid fc Spergell ()2009[ ) constrain the LRG 
HOD with the counts-in-cylinders multiplicity function 
and the correlation function through populating halos in 
a simulation. They obtain a good fit by using SO halos 
and the high-mass end slope is found t o be close to unity . 
Compared to friends-of-friends (FoF: [Davis et al.lll985f) 
finder, the SO halo finder does not have the problem of 
linking two halos by a thin bridge of particles. We plan 
to pursue analytic models of galaxy c lustering based on 
SO halo properties (Ti nker et al.l '2008) in future work. 

Because of the steep high mass slope, our model 
fits predict that massive clusters should host multiple 
LRGs. For example, the LRG occupation number of a 
2 X 10^^h~^M(7, cluster would be about ten. Such a pre- 
diction can be tested with a cluster catalog if the mass 
can be determined. Using a sa mple of X-ray s elected 
galaxy clusters at 0.2 < z < 0.6, iHo et al.l (|2007| ) assess 
cluster membership for LRGs based on their photometric 
redshifts and assign halo mass based on X-ray luminos- 
ity. They define halos as objects with mean density of 
200 times the critical density rather than the mean back- 
ground density as we do, and they assume flm = 0.238. 
After correcting the differences in the halo definition and 
cosmology, our HOD result for the faint LRG sample 
matches theirs in the regime of {N{AI)) ~ a few. (How- 
ever, near the cutoff mass they find a much larger oc- 
cupation number, while we have {N{M)) < 1.) In their 
catalog, there are only two very massive clusters with 
masses of ~ 8 x lO^^h'^MQ and - 1.1 x lO^'^h'^MQ 
(corrected to be consistent with our halo definition), and 
the (corrected) numbers of LRG members are about 8 
and 14, respectively. We also performed a rough calcula- 
tion to a ssociate LRGs in th e faint sample to MaxBCG 
clusters (Koes ter et al.l [20071 ) in the overlapped sky re- 
gion and redshift range. The mass of each cluster is 
estimated from the total number of MAIN galaxies in- 
side the virial radiu s, with calibration by weak lensing 
([Sheldon et al.|[2007l ). For each cluster, we infer the halo 
radius and velocity dispersion from the halo mass cor- 
rected to match our halo definition. LRGs that fall in 
the projected halo radius and three times the velocity 
dispersion along the line-of-sight direction are assigned 
as cluster members (with no completeness and edge cor- 
rection applied). The number of LRG members for clus- 
ters more massive than 6 x IQ^'^h^^ Mq is found to range 
from to 7. The predicted number of LRGs in high mass 
halos (Fig. [1]) appears approximately consistent with the 
estimates from clusters. However, given the uncertain- 
ties in the mass estimator and the small number statis- 
tics, more work is needed to test our derived values of 
the high mass slope. 

By modeling clustering of LRGs with different lu- 
minosities, we infer how the mean luminosity of cen- 
tral galaxies changes with the mass of their host halos, 
cx MP with p ~ 0.46-0.51 for M ~ lO^/i-iM©. 
This is consistent with the relation M (x L"^ found by 
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Fig. 6. — Modified projected two-point correlation function and six-parameter HOD fit. Panel (a) shows the origina l i»p(rp) data points 
(fill ed ci rcles) and the modified ones (open circles), which approximate those predicted by the semi-analytic model of IBower et al.l 1 120061) 
(seelXInicida ct al. 200^. The thick curve is the best fit to the original data from the five-parameter HOD model and the dotted curve is 
that from the HOD model presented in § 3.1. The thin curve is the best fit to the modified data from the six-parameter HOD model (see 
the text). The two curves are best HOD model fits. Panel (6) shows the best fit mean occupation functions. The thick curves are from 
fitting the original Wp{rp) with the five-parameter model and the thin ones are the Ax'^ < 1 envelopes of mean occupation functions from 
fitting the modified Wp{rp) with a six-parameter model (see the text) and and those for all galaxies are shaded. Dasficd, dotted, and solid 
curves are for central, satellite, and all galaxies, respectively. 



IMandelbaum et aTl (I2006f) based on galaxy-gala xy lens- 
ing measurements. iPadmanabhan et al.l (j2009( ) model 
the clustering of a sample of photometrically selected 
SDSS LRGs that are fainter than those in the spectro- 
scopic samples we use. Combining their HOD modeling 
results with ours for LRGs at z ~ 0.3, we find that the 
value of p is in the range of 0.4-0.5 over a larger range 
of LRG luminosity. From a study of the halo occupa- 
tion statistics of galaxies using galaxy groups identified in 
the Two Degree Fi el d Galaxy Redshift Survey (2dFGRS; 
IColless et al.][200l . lYang et al.l (l2005l) infer a L^-M re- 
lation that is well described by a broken power law with 
p ~ 2/3 and 1/4 below and above 10^^ h~^M(7), respec- 
tively. The inferred Lc-M relation by iVale fc Ostrikeil 
(|2004f ) from matching the luminosity function from the 
2dFGRS with the theoretical sub halo mass function has a 
similar high-mass slope p ~ 0.28. iMore et al.l (|2009f ) con- 
strain the high-mass slope to be 0.28^Q'Qg for z < 0.072 
SDSS galaxies based on satellite kinematics. The high- 
mass slope extrapolated from the HOD modelin g resu lts 
of SDSS MAIN galaxies (jZheng. Coil, fc Zehavill200l is 
also about 0.3. The value we obtain at the high-mass end 
from modeling the LRG clustering differs significantly 
from these results. It is interesting to see whether such 
a difference can be explained by the difference in the 
galaxy samples. The 2dFGRS luminosity is in the bj 
band, while we use the SDSS g-band luminosity. Since 
the wavelength coverages of these two bands are close, 
it seems unlikely that the band difference can cause the 
apparent discrepancy. The other thing to notice is that 
the 2dFGRS galaxies and the SDSS MAIN galaxies have 
a mean redshift ^0.1 and the LRG galaxies in our anal- 
ysis are located around redshift 0.3. Could the discrep- 
ancy indicate an evolution effect over the intervening ~2 
billion years? Through fitting restframe B-band galaxy 
luminosity functions at different redshifts using a CLF 



approach, iCooravl (|2005l ) finds that the data are compat- 
ible with a halo-mass-dependent central galaxy luminos- 
ity evolution, with the high m ass slope p incr easing with 
redshift. The fitting results of ICooravl ()2005D imply that 
p could be as high as 0.5 at z '--^ 0.3, close to our inferred 
value. 

If we take the inferred values of p at z ^ 0.3 from 
our analysis and at z ~ 0.1 from the 2dFGRS stud- 
ies at face value, they indicate that the luminosity evo- 
lution of central bright galaxies depends on halo mass 
in the sense that either galaxies in more massive halos 
fade more or the fraction of stars that were assembled 
into central galaxies from mergers betwe en z ^ 0.3 and 
z ^ 0. 1 is smaller in more massive halos. iBernardi et al.l 
(|2006D 's study of the properties of early- type galaxies in 
the SDSS as a function of local environment and red- 
shift suggests that star formation in early-type galaxies 
happens earlier in dense regions and lasts over a shorter 
time-scale, which implies that central galaxies in more 
massive halos on average experience star formation at 
an earlier epoch. Since younger stellar populations fade 
faster than older ones, this seems to rule out the possi- 
bility that luminosity of central galaxies in more massive 
halos fades more between z ~ 0.3 and z ~ 0.1. We 
are left with the possibility that mergers in more mas- 
sive halos are less efficient in adding stars to the central 
galaxies, which appears to be consistent with our finding 
of the drop of Mi /Mmin and its interpretation based on 
the competition between accretion and destruction as a 
function of halo ma ss. This is also in line with the results 
of LRG evolution in lBrown. Zheng. White, et al. I (j2008l ) 
from HOD modeling of their clustering from z ~ 0.2 to 
z - 1.0 in the NOAO Deep Wide-Field Survey. 

Since LRGs naturally separate out massive halos, our 
HOD modeling of the two-point cross-correlation func- 
tions between LRGs and i, galaxies circumvents the 
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usual challenge in constraining the HOD of galaxies 
in massive halos from their two-point auto-correlation 
function. We show that the cross-correlation data are 
consistent with the case where the distribution of L* 
galaxies inside massive halos roughly follows that of the 
matter at distances greater than ^ 0.2/i^^Mpc and that 
the mean number of L, galaxies scales with halo mass 
as M^-^. The slope of the mean occupation function is 
steeper than what is found for z ~ galaxies, which is 
^ 1.10. There may not be inconsistency between the 
results. The constraint for the slope for z ^ galax- 
ies from auto-correlation function is not sensitive to the 
occupation in very massive halos, while here from cross- 
correlation with LRGs, we have a better constraint on the 
occupation function in massive halos. On the other 
hand, the selection of galaxies and its redshift depen- 
dence make it hard to do a precise comparison between 
results at z ~ and at z ~ 0.3. 

The luminosity dependence of the LRG-L, galaxy 
cross correlation depends on scale in a rather complex 
way (jEisenstein et al.l l2005af ). By separating contribu- 
tions from pairs of galaxies with central and satellite 
LRGs in common halos and in different halos, our HOD 
modeling results explain these trends, in a relatively 
transparent way. At a fixed scale, the luminosity de- 
pendence in the cross correlation reflects the luminosity- 
dependent HOD of LRGs. As the scale in considera- 
tion changes, the relative contributions of the one-halo 
central-satellite, the one-halo satellite-satellite, and the 
two-halo LRG-L* pairs to the cross-correlation function 
vary, which leads to the scale-dependent luminosity de- 
pendence of the cross correlation. 

Our LRG modeling results establish the relation be- 
tween massive galaxies and dark matter halos at z ^ 0.3, 
which itself provides useful tests to models of formation 
of massive galaxies. lAlmeida et al.l (2008) present pre- 
dictions for properties of LRGs in semi-an alytic galaxy 
forma tion models. They show that the iBower et al.l 
(|2006f ) model, whic h is based on the Millennium simu- 
lation (jSpringel et al. 2005), predicts a z = 0.24 LRG 
luminosity function that is in good agreement with the 
observation, although it fails at z = 0.5 . Compared to 
the measurement, the lBower et al.l (|2006[ ) model seems to 
predict the z = 0.24 LRG two-point correlation function 
remarkably well on both small and large scales (see their 
Fig. 13). However, the LRG HOD predicted in this model 
differs significantly fro m our modeling res ults presented 
in this paper: (a) in the lBower et al.l (|2006l ) model, LRGs 
in the fainter sample can reside in halos of 10^^/i~^Mq 
(with mean occupation number of ~0.01), much lower 
than the mass scale we infer; (b) the mass of halos that 
can on average host one LRG is about 3 xl0^^h_^_M^ 
highe r than what we find (see Fig. 1 1 in lAlmeida et al.l 
[lOOland Fig.16 in lWake et all 120081 ): (c) the probabil- 
ity distribution of LRG host halo masses is broad (see 
Fig.16 in I Wake et al]l2008l ). ranging from lO^^/j-i^f^ 
10^^h~^MQ, rather than narrowly peaked around a few 
times lO^^h-^MQ as we find (Fig.[TJi). 

Does the discrepancy between the theoretically pre- 
dicted HOD and our observationally inferred HOD im- 
ply that our HOD parameterization is not generic enough 
to model LRG clustering? To investigate the problem, 
we modify the five-parameter HOD model (Appendix |BJ 



to mimic th e shape of t h e me an occupation function 
predicted in IBower et al.l (|2006f) model. The parame- 
terization for the satellite mean occupation function re- 
mains unchanged (with three parameters). For the cen- 
tral galaxy occupation function, we model it as a lin- 
ear (in logarithmic space) ramp going from (iVmin) at 
Maihi to unity at M„ and staying at unity for M > Mu- 
In total, this parameterization has six free parameters. 
With the measured Wp(rp), we find that the best- fit HOD 
from this six-parameter model (not shown in Figure [H)) 
closely follows the result from the five-parameter model 
(shown as thick curves in Figure [6l)). Therefore, change 
in the HOD parameterization does not solve the discrep- 
ancy and the five-parameter model is not inadequate in 
m odeling LRG c l usterin g. A close look at Figure 13 
in lAlmeida et al.l ([20081 ) shows that the predicted two- 
point correlation function does not match the data per- 
fectly — it is ~15% higher on scales of 0.1-2/i~^Mpc 
and ~ 30% lower on scales larger than 2/i~^Mpc (note 
that the vertical range of the plot is over eight orders 
of magnitude). We therefore modify the amplitude of 
the observed w.p{rp) data points to mimic the predicted 
correlation function and perform an HOD fit with the 
six-parameter model (by adopting the predicted number 
density, which is 10% lower than the observed one). The 
results on the mean occupation functions are shown in 
panel {b) of Figure[6]as thin curves. Note that the shaded 
region represents the Ax^ < 1 envelopes and the linear 
ramp for the central galaxy mean occupation number 
(thin dashed curves) has not reached unity at the high- 
est mass in the plot. The mean occupation function can 
extend to halos of mass as low as 2 x \Q^^h~^MQ and 
reaches unity around 4 x \Q^^h~^MQ, which appears to 
be approximately c onsistent with the prediction of the 
IBower et al.l ()2006[ ) model. From the above investiga- 
tions, we conclude that the discrepancy between the the- 
oretically predicted LRG HOD and the observationally 
inferred LRG HOD refiects the imperfection of the semi- 
analytic galaxy formation model — the 15 — 30% discrep- 
ancies with observed clustering are real and physically 
significant — rather than limitations in our HOD pa- 
rameterization. The results suggest that the mechanism 
of turning blue galaxies to red in the semi-analytic model 
is too efficient in halos of a few times 10^^/i~^Mo and 
not efficient enough in more massive halos. Our HOD 
modeling results thus provide important tests to galaxy 
formation theory. 

In combination with pas sive evolution of LRGs and a 
halo merger history (e.g., IWhite et al.l I2007t ISeo et all 
l200l . our modeling results can be used to predict the 
HOD of LRGs and the clustering of LRGs at lower or 
slightly higher redshifts. Supplemented with correspond- 
ing observations at these redshifts, we would be able to 
test our understanding of the formation and evolution 
of massive galaxies. Constraints on the HOD of LRGs 
are also useful for cosmological parameter investigations 
based on LRG clustering. For example, the most precise 
measurements of the large scale galaxy power spe ctrum 
have come from the SDS S LRG sample ( Teg mark et al.l 
120061 : [Percival et al.l2007f ). and the principal limitation in 
interpreting these measurements is the uncertain level of 
scale-dependent bias between galaxy and matter power 
spectra in the mildly non-linear regime. With HOD 
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constraints like those derived here, thi s scale-dependen t 
bias can be calculated and corrected (|Yoo et al.ll2009f ). 
As another example, HOD constraints could be com- 
bined with galaxy-galaxy lensing measurements of the 
LRG-mass cross-correlation (R. Mandelbaum et al., in 
preparation) to improve determinations of erg and flm 
lYoo et all 120061 ). The role of LRGs in cosmological 
studies seems destined to grow with surveys that target 
large numbers of LRGs to measure baryon acoustic_os- 
cillations, including AAOmega LRG survey (jRoss et alj 
[26O8) and the Baryon Oscillation Spectroscopic Survey 
(BOSS), part of a proposed successor to SDSS-II. Under- 
standing the evolving relation between LRGs and dark 
matter halos will be crucial to exploiting their power as 
cosmological probes and to revealing the physics that 
governs the formation of the most massive galaxies in 
the universe. 
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APPENDIX 

ROLE OF HALO MASS SCALES ON THE DEPARTURES FROM A POWER LAW IN THE GALAXY TWO-POINT 

CORRELATION FUNCTION 

Departures of the ga laxy two-point corre l ation function from a pure power law have been observed for both low and 
high redshift galaxies (Hawkins ct al. 26631: [Zehavi et al.ll2004l : lOuchi et al.ll2005l : iCoil et al.ll2006HLee eraI1l2006f l. In 
the two luminosity-threshold LRG samples studied in this paper, t he departures are al so clearly seen, shown as an 
upturn in the two-point correlation function at small scales (see also lZehavi et al.ll2005aD . Such departures have been 
nicely explained by the HOD model as the t ransition froru a reg ime dominated by one-halo pairs on small scales to that 
dominated by two-halo pairs on large scales (jZehavi et al.ll200"^ . The strength of the departures of the galaxy two-point 
correlation function from a pure power law is closely related to the amplitude of the one-halo term, which itself depends 
on the scatter in the occupation number and the halo mass function (Benson ct al. 2000; Bcrlind & Weinberg 2002). 
It would be helpful to gain a better understanding of the key factor that determines the strength of the departures. 

For luminosity-threshold samples, there are two mass scales, the characteristic minimum mass Mmin of halos that 
host central galaxies above the luminosity threshold and the mass scale Mi of halos that on average host one satellite 
galaxy above the luminosity threshold. An additional mass scale is the nonlinear mass M^i, marking a transition in 
the halo abundance from a power-law form to an exponential cutoff at high mass. We can thus identify two ratios 
that can shape the one-halo term: Afi/Af,„in and M^in/Mn\- The Afi/Mmin ratio affects the shape of the galaxy 
occupation functions and tells us how quickly the transition from sub-Poisson to Poisson scatter occurs when going to 
higher halo masses. A smaller Mi/Mj^in increases the importance of one-halo pairs by increasing the satellite fraction 
and thus results in a higher amplitude of the one-halo term. The Afmin/Mni ratio determines which part of the halo 
mass function the galaxy sample probes. The slope of the one-halo term reflects the drop of the halo mass function 
toward high masses. For a galaxy sample that probes the exponential tail of the halo mass function (Afmin/Mni > 1), 
we expect a steep drop in the one-halo term and thus a sharp upturn around the one-halo to two-halo transition scale. 

To figure out the relative importance of Mi/M^i^ and Mmin/Mni on the departures from a power law in the galaxy 
two point correlation function, we calculate the predicted correlation functions on a 3 x 3 grid of Mi/M^in and 
Mmin/Mni, as shown in Figure 1X1 In each panel, the short dashed curve is the one-halo term from central-satellite 
pairs, the long dashed curve is the one-halo term from satellite-satellite pairs, and the dot-dashed curve is the two-halo 
term. The dotted line is a power-law fit to Wp{rp) in the range of l-10/i~^Mpc for comparison. The satellite-satellite 
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Fig. 7. — Impact of the Mi/Mmin and M^in/M^i ratios on departures from a power law in the galaxy two-point correlation function. 
Each row (column) has the same Mi/Mmin (Mmin/A^ni) with the value marked on the right (top) of the plot. In each panel, the solid line 
shows the predicted projected correlation function, Wp(rp). The short and long dashed curves are the one-halo terms from central-satellite 
and satellite-satellite galaxy pairs, respectively, and the dot-dashed curve is the two-halo term. The dotted line is a power-law fit to Wp{rp) 
in the range of l-10/i~^Mpc to guide the eye. 

pair contributions are similar at a fixed Mi/M„iin, since (NgutiM)) is assumed to be proportional to halo mass. In 
a massive halo of fixed mass, however, the contribution of central-satellite pairs relative to satellite-satellite pairs is 
larger for a higher Mj^in/Mni sample because of its low satellite occupation number. This leads to a larger slope across 
the central-satellite and satellite-satellite contributions for a higher M^i^/M^i sample, thus a steeper upturn in the 
correlation function^. In more detail, the slope across the central-satellite and satellite-satellite contributions should 
depend on the Mi mass scale, but Mi/Mmi„ should not vary by an extremely large factor at least for a threshold 
galaxy sample. So a dependence on Mi can be translated to that on Mmin/Mni. When probing the exponential tail 
of the halo mass function for observed galaxy samples, we expect a steeper slope in both the central-satellite and 
satellite-satellite one-halo terms. This in turn leads to a more prominent rise at small scales. 

From the results, we conclude that overall the Afniin/Mni ratio plays a much more important role than Mi/Mmin 
in driving departures from a power law in galaxy two-point correlation functions. That is, a steeper drop in the halo 
mass function is translated to a steeper radial cutoff in the one-halo term, leading to a stronger inflection where the 
one-halo and two-h alo term join. This explains why the departure is stronger for more luminous galaxy samples (e.g., 
iZehavi et al.ll20d5ai rbl) that probe the high mass end of the halo mass function. This also explains why the departures 
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Fig. 8. — Dependence of the LRG HOD on cosmology. Panel (a) shows the measured two-point correlation function (data points with 
error bars) of the Mg < —21.2 sample and the best HOD fits (solid curves) for cosmological models differing only in erg. Note that the best 
fits overlap with each other. Panel (6) shows the mean occupation functions (solid curves) corresponding to the best fits, separated into 
central (dotted) and satellite (dashed) contributions. From left to right, the HODs correspond to eg values increasing from 0.65 to 1.00 
with a step-size 0.05. The two bottom panels are similar, but for the Mg < —21.8 LRG sample. 



are more prominent at higher redshifts (e.g., lOuchi et al.ll2005fl — the fast drop of the nonhnear mass Mni toward 
high redshift makes the observed (bright) galaxies more hkely to probe the exponential tail of the halo mass function. 
A more detailed theoretical investigation of the departure from a power law in the two-point correlation function of 
galaxies and its dependence on galaxy properties will be presented elsewhere. 

LRG HODS WITH A FIVE-PARAMETER MODEL AND DEPENDENCE ON COSMOLOGY 

The LRG clustering modeling in the main text assumes a flexible HOD parameterization, which reveals details of the 
constraining power of the two-point correlation functions on the HOD. In general, the HODs for luminosity-threshold 
samples predicted by galaxy formation models can be well described by a less flexible (i.e., more restricted) parametric 
form. The mean occupation functi on is well charac terized by a step-like fu nction for central galaxies and a power law- 
like function for satellite galaxies (jKravtsov et al.l 2004; Z heng et al.ll2005h . We present here the modeling results for 
the two luminosity-threshold LRG samples with an alternative, five-parameter HOD model and give their dependence 
on cosmology. This set of results would be useful for comparisons with other work and for making mock catalogs in a 
wide range of cosmological models. 

The mean occupation function of a luminosity-threshold LRG sample, being the sum of central and satellite mean 
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occupation functions, is parameterized as (see IZheng et al.l[200l IZheng. Coil, k Zehavil[200l 



(iV(M)) 



1 + erf 



log M - log M„ 



ClogM 



M-Mo 
Ml 



where erf is the error function 



erf (x) 



dt. 



(Bl) 



(B2) 



The distribution of the occupation number of central galaxies and satellite galaxies are assumed to follow the nearest- 
integer and Poisson distributions, respectively, as usual. The five free parameters are the mass scale Mmin and 
width CTiogM of the cutoff profile for the mean occupation function of central galaxies and the cutoff mass scale Mq, 
normalization M[, and high mass slope a of the mean occupation function of satellite galaxies. 

We vary the normalization erg (at z = 0) of the matter fluctuation power spectrum from 0.65 to 1.00 with a step-size 
0.05. Other cosmological parameters are assumed to be {flm, ^A, ^b, ris, /i)=(0.24, 0.76, 0.04, 0.95, 0.73). We perform 
an MCMC run for each LRG sample under each cosmological model and obtain the marginalized distribution for each 
of the five HOD parameters. 

In general, for a higher erg, the nonlinear mass increases and there are more high mass halos. To conserve the galaxy 
number density, the cutoff mass scale Afmin needs to increase. The increase in Mmin turns out to be slower than that 
in Af„i, which leads to a lower halo bias that is necessary to maintain the large scale clustering of galaxies (i.e., the 
square of bg(T%). The width triogM of the cutoff increases so that some LRGs are populated to lower mass halos to 
further adjust the large scale bias factor and the galaxy number density. The mean occupation function of satellites 
also responds to the (j% change to match the small scale clustering amplitude: the mass scale M[ increases and the 
slope a decreases with increasing erg. For a large range o f gg , our HOD mod e ling y ields almost identical best fits to 
the data points, similar to what is found and discussed in IZheng fc Weinbergl (|2007l ). 

While the dependence of log Mmin, ciogM, logMo, logM{ or a on erg appears to be quite close to hnear, we fit the 
results by adding a quadratic term for a better accuracy. The HOD parameters for the Mg < —21.2 LRG sample are 

log A/min = 13.673 -f 1.419(crg - 0.8) - 1.706(crg - 0.8)^ 
(TiogM = 0.621 + 0.908(crg - 0.8) - 0.935(crg - 0.8)^ 
log Mo = 12.339 + 0.658(crg - 0.8) + 9.206(crg - 0.8)^ 
logM{ = 14.533+ 1.248(crg - 0.8) - 1.394(crg - 0.8)^ 
a= 1.832 - 1.326(crg - 0.8) + 0.523(crg - 0.8)^. 



(B3) 



The typical 1-cr uncertainties for the five HOD parameters are 0.06, 0.07, 0.6, 0.025, and 0.08, respectively. The HOD 



parameters for the Mg < 



-21.8 LRG sample are 

log Mmin = 14.304 - 
criogM= 0.797- 
logMo = 12.491 - 
log M{ = 14.946 
a= 1.717 



1.694(crg - 0.. 

0.761(crg-0.. 
1.476(cr8-0.. 
1.616(cr8-0.. 

0.589(crg - 0.. 



- 1.810(crg - 0.8)' 

- 0.614(crg - 0.8)' 
29.983(crg - 0.8)' 

- 0.712(crg - 0.8)' 

- 7.437(crg - 0.8)' 



(B4) 

0.1, and 0.4, respectively. All 



The typical 1-cr uncertainties for the these five HOD parameters are 0.06, 0.055, 0.7 
the masses are in units of h~^MQ. 

For results with cosmolog ical models with 0^ different than 0.24, on e only needs to change the three mass scales 
by a factor of f7,„/0.24 (see IZheng eralll200l IZheng k Weinbe rg"2007). This scaling assumes that the shape of the 
linear power spectrum stays fixed, with a change in h or compensating the impact of changing flm- Therefore, the 
solutions here cover a wide range of cosmological models in the (f2„i,crg) plane. 

HALO NUMBER FLUCTUATION IN THE LRG SURVEY VOLUME AND THE COVARIANCE MATRIX 

In our HOD mo deling of the LRG two-point correlation functions, we use a theoretical halo mass function 
penkins et al.|[200li) . which is a fitting formula based on A^-body simulations. In the LRG survey volume, which is 
0.72/i^Gpc~^ (|Eisenstein et al.ll2005bl ). fiuctuations in the number of massive halos (e.g., with mass above 10^^h~^ Mq) 
are expected. Since LRGs reside in massive halos, the details of their small scale clustering may be sensitive to such 
fluctuations. In this Appendix, we investigate whether these fluctuations introduce any systematic effect on the HOD 
modeling. 

iJing et al] (|2007f ) performed A^-body simulations with 1024^ particles in a box of 1800/i^^Mpc on a side. For our 
investigation, we make use of the catalog of massive halos (> 10^^h~^MQ) identified in the z = 0.274 outputs from 
four independent simulations. The volume of each realization is divided into eight octants. Each octant has a volume 
similar to the LRG survey volume, so altogether we have 32 sub-volumes to investigate the mass function variations. 
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Fig. 9. — (Left): Halo mass function and its scatter from the IJine et ai.l II2007I1 simulation, for a volume similar in size to the LRG 
survey volume. The green solid curve is the mean halo mass function in the simulation, the green dashed curves mark the scatter and the 
red dashed curves denote a Poisson scatter around the mean. The black solid, dotted, and dashed curves are the Jenkins, Sheth-Tormen, 
and Press-Schechter mass functions, respectively (see the text). (Right): Contribution of LRGs in massive halos to small-scale Wp and the 
variation due to the fluctuation of halo mass function in a finite volume. The green solid curve is the best fit from a five-parameter HOD 
model to the measurements (green points). The green dotted and dashed curves are the one- halo and two-halo terms, respectively. The 
red solid curve is Wp calculated by excluding halos more massive than lO^^h-'^MQ. The dotted red curve shows the one-halo term from 
halos below 10^^h~^MQ, and the blue curve shows it for halos above this value. The black solid curves are the one-halo term contributions 
from halos more massive than 10^^h~^MQ in each of the 32 simulation subvolumes that have the same size as the LRG survey volume. 
The bottom panel shows the comparison of the diagonal jackknife errors and uncertainties introduced by the fluctuation of the number of 
massive halos in the LRG survey volume. See text for more details. 

The left panel of Figure [C] shows the fluctuation of the number of massive halos among the 32 sub- volumes. The 
green solid curve is the mean halo mass function in a sub- volume and the green dashed curves mark the scatter around 
the mean. The scatter turns out to closely follow that of a Poisson distribution (the two dashed red curves). The 
fluctuation in the number of halos is about 15% for M > 10^^ Mq and increases with halo mass. For comparison, 
we also plot three frequentl y used theoretical f unctions for the same cosmological model assumed in the simulation. 
The Jenkins mass function (| Jenkins et al.|[200ll ). shown as the black solid curve, gives a good description of the mean 
halo mass function at 10^^ M p, and becomes a factor of two higher at 2 x 10^^/i~^Mq. The Sheth-Tormen mass 
function (Shoth fc Tormenlll999t black dotted curv e) is slightly higher than the Jenkins mass function, while the Press- 
Schechter mass function ( Press fc Schechterl 119741 : black dashed curve) underestimates the mass function by a factor 
of five in the mass range considered here. 

The right panel of Figure [C] shows the effect of the fluctuation in the number of massive halos on the small scale 
clustering of the Mg < —21.2 LRGs. The green solid curve is the best fit to the measured Wp (green points) from 
the five-parameter HOD model (Appendix |B]) with the same cosmological model as used in the simulation. The green 
dotted and dashed curves are the contributions from the one-halo and two-halo term, respectively. With the best-fit 
HOD model fixed, the red solid curve show the predicted Wp when only keeping halos less massive than 10^^h~-^ Mq 
in the calculation. The red dotted curve is the one-halo term from halos less massive than 10^^h~^MQ. 

The blue solid curve is the mean contribution to the one-halo term from halos of M > 10^^h~^ Mq. The fluctuation 
in the number of massive halos would lead to a fluctuation around this mean contribution. To see this, we populate the 
massive halos in the 32 sub-volumes according to the best-fit HOD model and measure the one-halo Wp in each sub- 
volume. Black solid curves show the measurements for individual sub- volumes. Such a fluctuation in the one-halo term 
introduces uncertainties in the small scale clustering. By adding the one-halo term from halos with M < 10^^h~^MQ 
(red dotted), that from halos above 10^^h~^MQ in each sub- volume (black solid) and the two halo-term (green dashed), 
we obtain the black points with error bars reflecting the fluctuation in the number of massive halos. On small scales, 
these error bars appear to have similar (somewhat smaller) amplitude to those from the diagonal elements of the 
covariance matrix (see the comparison in the lower panel), which is computed through jackknife technique. 

We see that the fluctuation of the number of massive halos in the LRG survey volume can be large, e.g., ~15% for 
M > 10^^/i~^Mo halos and ~50% for M > 2 x IO^^H'^Mq halos. This introduces a fluctuation in the smafl scale 
clustering of LRGs. However, our investigation suggests that such a fluctuation in Wp should be correctly reflected in 
the jackknife covariance matrix. Therefore, it is sufficient to use the mean mass function in modeling the two-point 
correlation functions of LRGs. 
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